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Abstract 

Graphene has demonstrated great promise for future electronics technology as well as fundamen- 
tal physics applications because of its linear energy-momentum dispersion relations which cross at 



the Dirac point 



However, accessing the physics of the low density region at the Dirac point has 



been difficult because of the presence of disorder which leaves the graphene with local microscopic 
electron and hole puddles 3" 5], resulting in a finite density of carriers even at the charge neutrality 
point. Efforts have been made to reduce the disorder by suspending graphene, leading to fabri- 
cation challenges md «c,e devices w hidl mak e .oca, „c P ic _e, rtS difficult flfl. 
Recently, it has been shown that placing graphene on hexagonal boron nitride (hBN) yields im- 



proved device performance 



8(. In this letter, we use scanning tunneling microscopy to show that 



graphene conforms to hBN, as evidenced by the 
images. However, contrary to recent predictions [9j, 



presence of Moire patterns in the topographic 



ipj , this conformation does not lead to a sizable 



band gap due to the misalignment of the lattices. Moreover, local spectroscopy measurements 
demonstrate that the electron-hole charge fluctuations are reduced by two orders of magnitude 
as compared to those on silicon oxide. This leads to charge fluctuations which are as small as 



in suspended graphene 



opening up Dirac point physics to more diverse experiments than are 



possible on freestanding devices. 
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Graphene was first isolated on si 
regions using an optical microscope 



icon dioxide because of the ability to image monolayer 



111 ] . However, the electronic properties of SiO"2 are not 



ideal for graphene because of its high roughness and trapped charges in the oxide. These 
impurity-induced charge traps tend to cause the graphene to electronically break up into 
electron and hole doped regions at low charge densit y w hich both limit device performance 
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131 ] . In order to create devices with less 



and make the Dirac point physics inaccessible {3- 5, 
puddles, the substrate must be removed or changed. One possibility to get rid of substrate 
interactions is to suspend graphenejf], Q], as shown by the drastic improvement in mobil- 
ity which has enabled the observation of the fractional quantum Hall effect in suspended 
devices [bll. However, the freely suspended monolayers are very delicate, leading to fab- 



rication difficulties as well as strain 



161 ] . Because of these difficulties, there have been no 



STM spectroscopy measurements of the local electronic properties of suspended graphene 
devices. All of this points to the need for new substrates that offer mechanical support 
to the graphene without interfering with its electrical properties. Recently, such a candi- 
date substrate has been found with the demonstration of high-quality graphene devices on 
hexagonal boron nitride (hBN)js|. Hexagonal boron nitride has the same atomic structure 
as graphene, but with a 1.8% longer lattice constant 17|, and shares many similar properties 
with graphene except that it is a wide-bandgap electric insulator 18J . The planar structure 
of hBN cleaves into an ultra-flat surface and the ionic bonding of hBN should leave it free of 
dangling bonds and charge traps at the surface resulting in less induced electron- hole pud- 
dles in graphene. Indeed, graphene on hBN devices exhibit the highest mobility reported on 
any substrate, as well as narrow Dirac peak resistance widths, indicating reduced disorder 
and charge inhomogeneity 8J . 

To study how the local electronic structure of graphene is affected by the hBN substrate, 
we prepare graphene on hBN devices for STM measurements. A schematic of the measure- 
ment set-up showing the graphene flake on hBN with gold electrodes for electrical contact 
is shown in Fig. [Tk. A typical STM image of the monolayer graphene showing the surface 
corrugations due to the underlying hBN substrate is shown in Fig. [lb. This image can be 
compared with an STM image of monolayer graphene prepared in a similar manner but on 
SiO"2 as shown in Fig. [lb. It is clear from these two images that the surface corrugations 
are much larger for graphene on SiO"2 as compared to hBN. This is due to the graphene 
conforming to the substrate and the planar nature of hBN as compared to the amorphous 
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SiC>2. Figure [TJd shows a histogram of the heights in the two images. In both cases, the 
heights are well described by gaussian distributions with standard deviations of 224.5 ± 0.9 



pm for graphene on SiC>2 and 30.2 ± 0.2 pm for grap 



on Si02 are similar to previously reported values 
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on hBN is similar to graphene on mica or HOPG2l|. Reducing the surface roughness is 



rene on hBN. The values for graphene 



20] while the distribution for graphene 



critical for graphene devices because local curvature can lead to electronic effects such as 



doping 



22] and random effective magnetic fields 23j . As the height variation of the graphene 



on hBN is as flat as HOPG, it has reached its ultimate limit of flatness. 
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Figure 1: Schematic device setup and topography comparison of graphene on hBN and Si02- 
(a) Optical microscope image of the mechanically exfoliated monolayer graphene flake with gold 
electrodes. The wiring of the STM tip and back gate voltage is indicated, (b) STM topographic 
image of monolayer graphene on hBN showing the underlying surface corrugations. The image is 
100 nm x 100 nm. The imaging parameters are tip voltage Vt = -0.3 V, tunneling current It = 100 
pA. (c) STM topographic image of monolayer graphene on Si02 showing significantly increased 
corrugations. The imaging parameters are tip voltage Vt = -0.5 V, tunneling current It = 50 pA. 
(d) Histogram of the height distributions for graphene on SiC>2 (blue squares) and graphene on 
hBN (red triangles) along with gaussian fits. 



Looking more closely at the topography of the graphene on hBN, we resolve its atomic 
lattice and also observe longer periodic modulations which result in a distinct Moire pattern 
as seen in Figs. [2^ and|2fc- These two images were acquired from different areas of the same 



graphene flake and show Moire patterns of different length due to a shift in the alignment 
of the graphene with the underlying hBN lattice. In the case of Figs. [2b and [2b, we find a 
long wavelength modulation of 2.6 nm and 1.3 nm respectively. By examining the Fourier 
transform of these images, we can learn about the underlying hBN substrate as well as its 
relative orientation with respect to the graphene lattice. In both cases, we see two distinct 
sets of peaks in the FFT, there is a set of six points near the center of the images which 
correspond to the Moire pattern and there are six additional points near the edge of the 
images which correspond to the graphene atomic lattice. The first observation is that the 
locations of these points is rotated with respect to each other in the two images. For Fig. 
[2b, the atomic lattice is rotated by 12.6 ± 1.0° from the horizontal. In the case of Fig. [2H, 
the atomic lattice is rotated by 18.5 ± 0.6°. In contrast the Moire pattern is rotated by 
29.9 ± 0.1° from the horizontal in Fig. [2b and 38.4 ± 0.2° from the horizontal in Fig. (2U. 
From the lengths of the Moire patterns and their angles, we can calculate the orientation of 
the graphene lattice with respect to the underlying hBN (see Supplementary Information). 
In the case of Fig. [2h, we find that the hBN is rotated by -5.4° from the graphene. On 
the other hand, for Fig. [2b, we find that the hBN is rotated by -10.9°. This difference of 
5.5° matches the difference in the orientation of the graphene lattices in the two images. 
Therefore, we conclude that the underlying hBN substrate is continuous and the graphene 
above it sits at two different angles. Atomic force microscopy images show that graphene on 
hBN tends to form flat regions separated by ridges and pyramids. As the two STM images 
were taken from different sides of one of these ridges, it is clear that the graphene can change 
orientation across these ridges. 

The scanning tunneling microscope is not only able to acquire images of atoms but can 
also map the local density of states. We have performed scanning tunneling spectroscopy 
of the graphene on hBN. Figure [3b, shows a typical dl/dV spectroscopy curve which is 
proportional to the local density of states. The curve is nearly linear in energy with a 
minimum near zero tip voltage indicating the energy of the Dirac point. The location of this 
minimum can be varied by applying a back gate voltage, V g , to the sample. The voltage on 
the back gate induces a charge on the graphene of n — aV g with a = 7.2 x 10 10 e/cm 2 V 
based on a parallel plate capacitor model. In the model, we have taken 285 nm of SiC>2 with 



a dielectric constant of 3.9 and 14 nm of hBN with a dielectric constant of 3-4 



24). There is 



a small uncertainty in the value of a of about 1% due to the unknown precise value of the 
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Figure 2: Real space and Fourier transforms of Moire patterns (a) STM topography images of a 
Moire pattern produced by graphene on hBN. The scale bar is 2 nm. The inset is a zoom in of a 
2 nm region and the scale bar is 0.3 nm. The imaging parameters are Vj = -0.3 V and It = 100 
pA. (b) Fourier transform of (a) showing the six graphene lattice points near the edge of the image 
and the long wavelength Moire pattern near the center of the image and around each lattice point. 
The scale bar is 10 nm -1 . The inset is a zoom in around one of the lattice points. Its scale bar 
is 2 nm -1 . (c) STM topography image from another region of the same graphene flake showing 
a different Moire pattern. The scale bar is 2 nm. The inset is a zoom in of a 2 nm region and 
the scale bar is 0.3 nm. The imaging parameters are Vt = -0.3 V and It = 100 pA. (d) Fourier 
transform of (c) showing the atomic lattice as well as the Moire pattern. The scale bar is 10 nm -1 . 
The inset is a zoom of the Moire pattern and its scale bar is 4 nm -1 . 

hBN dielectric constant. Figure [3b plots dl/dV as a function of tip voltage and gate voltage. 
The white line follows the minimum in the dl/dV curves for each value of the gate voltage. 
This minimum occurs when the Fermi energy of the tip lines up with the Dirac point. We 
observe that the location of the minimum changes more quickly when the Dirac point is 
near zero tip voltage which is consistent with the linear band structure of graphene. There 
is also a dark ridge that occurs at decreasing tip voltage as the gate voltage is increased. 
This is due to the effect of the voltage on the tip acting as a local gate and changing the 
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density of electrons in the graphene 25 ] . 
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Figure 3: Spectroscopy of graphene on hBN as a function of gate voltage (a) dl/dV spectroscopy 
showing a nearly linear density of states as a function of energy (tip voltage), (b) dl/dV spec- 
troscopy as a function of tip voltage and gate voltage. The white line corresponds to the minimum 
in the dl/dV curves and represents the Dirac point, (c) Energy of the Dirac point as a function of 
gate voltage. The red curve is a fit assuming a linear band structure, (d) Energy versus momentum 
dispersion relations for the case of graphene and hBN having the same lattice constant and zero 
angle mismatch (black curve) and two curves with 1.8% lattice mismatch. The blue curve has 
-5.45° angle mismatch and the red curve has -10.9°. 



Our local spectroscopy measurements indicate that there is no band gap induced in 
graphene on hBN, not even locally. These results disagree with earlier theoretical calcula- 
tions which predicted the opening of a band gap of order 50 meV when graphene is placed 
on hBN, because of the breaking of sublattice symmetry (9), [loj]. This discrepancy is ex- 
plained by the 1.8% mismatch between the graphene and hBN lattices and the different 
orientations of the two lattices, which were both neglected in Refs. {^[loj]. Taking them into 
account, one expects that, while one of the carbon atoms may sit over a boron (nitrogen) 
atom at one location, this alignment gets lost a few lattice constants away. In large enough 
systems, carbon atoms should therefore have the same probability to have a boron or a 
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nitrogen atom as nearest neighbor in the hBN layer, regardless of their sublattice index. 
We numerically checked the validity of this hypothesis by calculating the interlayer hopping 
potential ~ 7j_ exp[— |r — r'|/£] from a carbon atom at r in the graphene layer to a boron or 
nitrogen atom at r' in a rotated hBN layer. We restricted ourselves to nearest- and next- 
nearest-neighbor interlayer hopping and chose the parameters 7 1 = 0.39 eV and £ = 0.032 
nm to fit known values for these hoppings in graphene bilayers 261 ] . We found that the 1.8% 
lattice mismatch alone is sufficient to make the hopping strength from a carbon atom to a 
boron or a nitrogen atom independent of the graphene sublattice index for systems of a few 
hundred unit cells, going down to a few tens of unit cells when the lattices are misaligned 
by about one degree. We incorporated the Fourier transform of this hopping potential into 
the low-energy Hamiltonian for graphene on hBN to find the energy-momentum dispersion. 
The inter-layer coupling is nonzero only for k = as well as for six additional vectors k 
associated with the Moire pattern. Most importantly, we found that the coupling between 
the A and B atoms in the graphene lattice with the boron and nitrogen atoms in the hBN 
are almost identical. Thus sublattice symmetry is restored and a gapless Dirac spectrum 
is recovered, albeit at slightly shifted values of K. This is illustrated in Fig. |3H. While 7_i_ 
depends in principle on whether hopping occurs between a carbon and a boron or nitrogen 
atom, we note that this does not break sublattice symmetry, and we checked that the spec- 
trum remains gapless, even when this discrepancy is taken into account. More details of our 
numerical approach are given in the Supplementary Information. 

By determining the energy of the Dirac point as a function of gate voltage, we can measure 
the Fermi velocity of electrons and holes in graphene. Figure [3b shows the energy of the 
Dirac point as a function of gate voltage. Graphene has a linear dispersion relation such 
that E = hvpk where vf is the Fermi velocity. Since it is a two-dimensional material, the 
density of electrons is given by n = g s g v 7rk 2 /(2tt) 2 where g s and g v are the spin and valley 
degeneracy which are both 2 for graphene. Therefore, the Dirac point should depend on 



gate voltage as E = hvp ^naV g . The red curve is a fit to the data from which we can 
extract the Fermi velocity. We find that vp = 1.16±0.01 x 10 6 m/s for the electrons and 
vp = 0.94 ± 0.02 x 10 6 m/s for the holes. Moreover, we observe an asymmetry between 
the Fermi velocity for electrons and holes of about 25% depending on the Moire pattern 
observed. The shorter Moire pattern has a higher Fermi velocity for holes while the longer 
one has a higher Fermi velocity for electrons. The origin of this asymmetry is unclear but 
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may arise due to next-nearest neighbor coupling which are not taken into account in our 
model. 

One of the main advantages to using hBN as a substrate for graphene as compared to 
SiC>2 is the improvement in electronic properties of the graphene which is believed to be due 
to the lack of charge traps on the hBN surface. Figure shows the topography of graphene 
on hBN over a range of 100 nm. Note that the height variation is less than 0.1 nm over the 
range of the image as compared to typical values of nearly 1 nm for graphene on SiC>2. We 
have performed dl/dV measurements at 1 nm intervals over the entire area of Fig. [4^. For 
each of these dl/dV curves, we have found tip voltage of the minimum which corresponds 
to the Dirac point. The results are plotted in Fig. 0b. We have done a similar analysis for 
a 100 nm area of graphene on SiC>2 and the results are plotted in Fig. [4b. The red and 
blue regions correspond to electron and hole puddles respectively. It is clear from these two 
images that the variation in the energy of the Dirac point is much smaller on hBN. The 
spatial extent of each puddle is also much smaller in the graphene on SiC>2 consistent with 
an increased density of impurities 13| . 

We can further quantify the disorder in the graphene by looking at a histogram of the 
energy of the Dirac point, Fig. 0H. The main part of the histogram for the Dirac point 
energy on hBN is well fit by a gaussian distribution (red line) with a standard deviation of 
5.4 ± 0.1 meV. In addition, there is a small extra bump in the distribution from the hole 
doped region near the bottom right of Fig. Hb- In comparison, the distribution on Si02 is 
much broader with a standard deviation of 55.6 ± 0.7 meV. These distributions in energy 
can be converted to charge fluctuations using n = E^/ir(hvF) 2 - We find that the charge 
fluctuations in graphene on hBN are a n = 2.50 ± 0.13 x 10 9 cm -2 while they are more than 
100 times larger for graphene on SiC>2, cr n = 2.64 ± 0.07 x 10 11 cm -2 . Our measurements 



for the charge fluctuations on Si02 are consistent with previous single electron transistor 
and STM0, 5] measurements which established the presence of electron and hole puddles in 
graphene on Si02- Furthermore, our measurements for the charge fluctuations in graphene 
on hBN show a very similar value to values extracted from electrical transport measurements 
in suspended graphene samples j^] implying that using hBN as a substrate provides a similar 
benefit to suspending graphene without the associated fabrication challenges and limitations. 

We have demonstrated that graphene on hBN provides an extremely flat surface that has 
significantly reduced electron-hole puddles as compared to Si0 2 . By reducing the charge 
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Figure 4: Spatial maps of the density of states of graphene on hBN and SiC>2 (a) Topography of 
graphene on hBN. (b) Tip voltage at the Dirac point as a function of position for graphene on 
hBN. (c) Tip voltage at the Dirac point as a function of position for graphene on Si02- The color 
scale is the same for (b) and (c). (d) Histogram of the energies of the Dirac point from (b) as well 
as a gaussian fit. The inset shows the same data but also includes the histogram for Si02 shown 
in red. The scale bar in all images is 10 nm. 

fluctuations, the low density regime and the Dirac point can be more readily accessed. 
Moreover, hBN allows this low-density regime to be reached in a substrate-supported system 
which will allow atomic resolution local probes studies of the Dirac point physics. 

Methods 



Thin and flat few layer hBN flakes were prepared by mechanical exfoliation of hBN sing 



crystals on SiC^/Si substrates. The hBN growth method has been previously described 271 ] 



Exfoliated graphene flakes were then transferred to the hBN using Poly(methyl methacry- 
late) (PMMA) as a carrier Then Cr/Au electrodes were deposited using standard elec- 
tron beam lithography. The lithography process leaves some PMMA resist on the surface of 
graphene which is cleaned by annealing in argon and hydrogen at 350° C for 3 hours 28]. 
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The device was then immediately transferred to the STM (Omicron low temperature STM 
operating at T = 4.5 K in ultrahigh vacuum (p < 1CT 11 mbar)). Electrochemically etched 
tungsten tips were used for imaging and spectroscopy. All of the tips used were first checked 
on an Au surface to ensure that their density of states was constant. 

The dl/dV spectroscopy was acquired by turning off the feedback loop and holding the 
tip a fixed distance above the surface. A small ac modulation of 5 mV at 563 Hz was 
applied to the tip voltage and the corresponding change in current was measured using lock- 
in detection. We also measured dl/dV curves with 0.5 mV excitation and observed the same 
results. 
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Supplementary Information 



I. MOIRE PATTERNS 

A Moire pattern occurs when the atoms in the graphene layer form a super-lattice struc- 
ture with the atoms in the hBN layer. In this section, we derive the conditions which lead 
to a Moire pattern and therefore predict the angles and lengths of the Moire pattern. 

We consider two superimposed hexagonal lattices defined by 

r ilM (m, n) = na i+ + ma,_ + (/x - i)d, (1) 

with the layer index % — 1,2 for the graphene and hBN lattices, respectively, and the 
sublattice index \i — 1 (A sublattice), 2 (B sublattice). In the hBN layer, /i = 1 for nitrogen 
and 2 for boron. We define d = (0, — l/v^)^ as the vector connecting the two sublattices, 
with lattice spacing a^. The vectors lattice vectors associated with the hBN, & 2 ±, are rotated 
by an angle counterclockwise with respect to the graphene lattice ai± = (±1/2, \/3/2)ai 
and are longer by a factor a 2 /ai = 1.018. Therefore the hBN lattice vectors are given by 
a 2 ± = (± cos(7r/3 ± 0), sin(7r/3 =F 4>)) a 2 

An arbitrary point in the graphene lattice can be represented by two integers (n, m) 
which correspond to the number of lattice vectors in the a 1+ and ai_ directions respectively. 
Similarly, a position in the hBN lattice is represented by the integers (r, s). If the graphene 
and hBN lattice are AB stacked at a given position, they will be AB stacked again when 
nai + + mai_ = ra 2+ + sa 2 _. In terms of xy-coordinates this gives the conditions 

n — m = 2 — (r cos(7r/3 + 0) — s cos(7r/3 — 0)) 
ai 

V3(n + m) = 2— (r sin(7r/3 — 0) + s sin(7r/3 + 0)) 

A Moire pattern occurs when these equations can be satisfied for integer values of (n, m) 
and (r, s). For a given ratio of a 2 /ai these conditions will only hold for certain special 
angles. In terms of the values of (n,m), the length of the Moire pattern can be written as 
L = ai\/n 2 + m? + nm and it is at an angle 9 = tan~ 1 ( 2 ^^ t ) with respect to a 1+ . However, 
a near commensurate condition can always be found leading to a Moire pattern over some 
finite length. 

We have numerically created lattices to reproduce the images in Figure 2. The results are 
shown in Fig. [ST] The first set of images correspond to the data shown in Fig. 2a and 2b 
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in the main text. It was created by using a hBN lattice with a 1.8% longer lattice constant 
and rotated by = —5.4° counterclockwise with respect to the graphene lattice. From the 
FFT of the lattice, Fig. [STb . we see that it matches the experimental figure, Fig. 2b very 
well. To create the shorter Moire pattern we must use a larger rotation angle. We find the 
best match occurs at = —10.9°. In general, the Moire pattern gets shorter as the angle 
increases. 





Figure SI: Simulated real space and Fourier transforms of Moire patterns (a) Simulated lattice 
showing a Moire pattern produced by graphene on hBN. The scale bar is 2 nm. (b) Fourier 
transform of (a) showing the six graphene lattice points near the edge of the image and the long 
wavelength Moire pattern near the center of the image and around each lattice point. The scale 
bar is 10 nm -1 . (c) Simulated lattice showing a shorter Moire pattern. The scale bar is 2 nm. (d) 
Fourier transform of (c) showing the atomic lattice as well as the Moire pattern. The scale bar is 
10 nm -1 
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II. THEORY CALCULATIONS 



We restrict the interlayer hopping potential to nearest-neighbor and next-nearest- 
neighbor hopping, and evaluate 



V^(m,n) = 7j_exp[-|r lAt (m,n) - r 2v (m' ,n')\/£\ 



(2) 



with (m',n') labelling the nearest and next-nearest neighbor sites to (m,n). Note that if 
one is a boron atom, the other one is a nitrogen atom. In this way, four interlayer couplings 
between the two sublattices are defined. The parameters 7j_ = 0.39 eV and £ = 0.032 
nm are calibrated to fit the interlayer couplings in bilayer graphene 26j . While 7j_ should 



in principle depend on the sublattice index /x in the hBN layer, this has no influence on 
our main conclusion, that the graphene spectrum is effectively gapless due to the mismatch 
between lattices and their different relative orientations, because it does not break sublattice 
symmetry in the graphene layer. 

We consider system sizes up to 600 x 600 unit cells, which is more than sufficient to extract 
the periodicity of the Moire patterns shown in Figs. 2a and 2c. Fig. IS2h shows a color plot of 
Vaa(^,^) for an interlayer rotation angle of <fi = —5.4° corresponding to the sample shown 
in Fig. 2a. The potential has a clear hexagonal periodicity, with a wavelength of about 
10 lattice sites, in quantitative agreement with Fig. 2a. From that hexagonal pattern, the 
Fourier transform V^k) of V^,(m, n) generically exhibits seven peaks, including a large one 
at k = 0, and six secondary peaks at wave vectors kj, % = 1, . . .6 corresponding to the 
Moire pattern. This is shown in Fig. IS2b . with hexagonally distributed secondary peaks at 
positions reflecting the central pattern of Fig. 2b. 

The hopping potential induces interlayer coupling which we introduce in the low-energy 
Hamiltonian. We keep only the seven dominant transitions and because the hBN bands 
are several eV's away from the Dirac points and V^,(kj^o)/^V(ko) < 0.1, we truncate the 
Hamiltonian to a 16 x 16 matrix 

6 



H(q,q') 



tf c (q)+ J]#BN(q + k*) 



i=0 



6, 



q,q 



2> 

,i=0 



C,BN 



'kj) + h.c. 



6, 



q+k;,q' ) 



(3) 
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Figure S2: Real space and Fourier transforms of the interlayer hopping (a) Vaa potential for an 
interlayer rotation of (j> = —5.4°. The scale bar is 2 nm. (b) Fourier transform of (a) showing a 
central peak at ko = and six additional peaks due to the Moire pattern. The scale bar is 0.8 



nm . 



with all H J s being 2x2 matrices, 

7o/o(q) 

#c(q) = , 4a 

7o/ *(q) o 

e B 7 i/ 1 (R^ 1 (0)q) 
#B N (q) = , (4b) 

7l A*(R^(0)q) e N 

^C,BN(ki) =1 _ , (4c) 

V BA (ki) V BB {\Li) j 

with /i(q) = 1 + 2 cos(q x ai/2) exp(—i\/3q y ai/2) and parameters 7 = 3.16 eV, 71 = 2.79 



eV, ee = 3.34 eV and eN = —1.4 eV[lO|. The low-energy graphene dispersion is obtained 
via diagonalization of %. We find that a gap exists in the spectrum only when V^A(kj) — 
V^B(kj) 7^ for at least one /x, thereby breaking sublattice symmetry in the graphene layer. 
At experimentally relevant rotation angles of a few degrees, we find that max[V r A1I/ (kj) — 
V^v(kj)]/V Mi ,(kj) < 1CT 4 for 150 x 150 systems, decreasing with size to less than 10~ 5 for 
600 x 600 systems. Consequently, we get an upper bound of A < 10~ 6 eV for the excitation 
gap A between the two graphene bands for the largest systems we investigated. The latter 
are still smaller than the experimentally observed domains so that such a small gap cannot 
be resolved. 

Large values for max[V^ tl/ (kj) — V r At /^/(k i )]/V^^(k i ) are obtained only for small systems or 
if the lattice mismatch is neglected. This explains the 50 meV gap reported in Refs. 
which we qualitatively reproduced (see Fig. 3d). Even without relative rotation of the two 
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layers, we found that the lattice mismatch is sufficient to effectively close this gap already 
for systems of sizes 300 x 300. We finally note that the expected different hoppings on the 
boron and nitrogen atoms induce differences in Vaa vs - Vab & s well as in Vba vs - Vbb- 
However, such differences do not break sublattice symmetry and therefore do not open a 
gap. 
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